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ABSTRACT 

We determine the nuclear spectral energy distribution (SED) from the radio to the gamma- 
ray band for the FR I radio galaxy NGC 6251, by using both data from the literature and 
analyzing HST and X-ray archival data. In the logv — log(i^F^) representation, the SED has 
two broad peaks, and it is remarkably similar to those of blazars. We show that the low-energy 
peak can be explained in terms of synchrotron radiation, while the high energy peak is most 
plausibly produced by inverse Compton scattering. This brings direct support to the FR I-BL 
Lacs unification model. We model the overall emission in the frame of a synchrotron self-Compton 
scenario, which well describes the SED of BL Lacs. The model parameters we obtain confirm 
also quantitatively the FR I-BL Lac unification model and imply a rather small viewing angle 
to this source (9 ^ 20°). NGC 6251 is the second radio galaxy, in addition to Centaurus A, for 
which a similar analysis has been performed. A significant improvement with respect to the case 
of Cen A is the absence of obscuration in NGC 6251, which strengthens the overall result. 

Subject headings: galaxies: active — galaxies: nuclei 



2 ESA fellow 
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1 Based on observations obtained at the Space Telescope 
Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Incorporated, un- 
der NASA contract NAS 5-26555. 



— galaxies: jets — galaxies: individual (NGC 6251) 
1. Introduction 

The basic scheme for unification of low-luminosity 
radio-loud active galactic nuclei (AGN) assume 
that radio galaxies of the Fanaroff-Riley class I 
(FR I, Fanaroff & Riley 1974) and BL Lac objects 
differ only because of a different viewing angle to 
the relativistic jet (Urry & Padovani 1995 for a re- 
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view). In BL Lacs, the jet is seen at small angles 
to the line-of-sight, and the non-thermal emission 
produced at its very innermost base, on the scale 
of ~ 10 15 — 10 16 cm, dominates the observed radia- 
tion. In FR I radio galaxies the jet points at larger 
angles, and the "BL Lac" component should be 
beamed away from the observer. 

Clearly, the most direct implication of the uni- 
fied models is that the extended properties (e.g. 
radio and optical emission, environment) of the 
beamed sources and their putative parent popu- 
lation should appear the same. Therefore, in the 
absence of any direct detection of the nuclear emis- 
sion in FR I, the usual way these models are tested 
involves the comparison of such extended compo- 
nents (e.g Urry et al. 2000, Antonucci & Ulvestad 
1985, Kollgaard et al. 1996). 

However, faint unresolved optical nuclei have 
been discovered in FR I radio galaxies through 
HST observations of a complete sample of 3CR 
sources (Chiaberge et al. 1999). The presence of a 
strong correlation between these optical nuclei and 
the radio cores argues for a common non-thermal 
synchrotron origin of both components. There- 
fore, this discovery provides qualitative support 
for the low-luminosity radio loud AGN unification 
models. 

A more quantitative analysis showed that cores 
of radio galaxies are overluminous with respect to 
what is expected from mis-oriented BL Lac jets, 
for typical values of the bulk Lorentz factor (Chi- 
aberge et al. 2000). It is therefore plausible that 
the optical cores of radio galaxies and the radia- 
tion observed in BL Lacs are not produced in the 
same region of the relativistic jets: in BL Lacs, a 
"fast jet spine" with Lorentz factors of <~ 15 — 20 
dominates, while in FR I we might be observing 
emission from a slower (but still relativistic) shear 
layer. 

Modeling of the overall nuclear SED could pro- 
vide an even more stringent test for the unified 
scheme, as it represents a fundamental tool for 
obtaining information on the physical conditions 
of the emitting region, as it is usually done for 
blazars (e.g. Ghisellini et al. 1998). Unfortu- 
nately, the nuclear SED of radio-galaxies arc in 
general not sufficiently well sampled to follow di- 
rectly this approach and radio galaxies are gener- 
ally not detected in the gamma-ray band. 



The first radio galaxy (and, to the best of our 
knowledge, the only case to date) for which a de- 
tailed modeling of the overall nuclear SED, from 
the radio to the gamma ray-band, has been un- 
dertaken is Centaurus A, the closest radio galaxy 
(Chiaberge et al. 2001). The result of such work 
is that, when considered in the comoving frame of 
the source, the nuclear SED of Cen A can be re- 
produced adopting the same emission model and 
physical parameters typical of BL Lacs of the same 
total power, strongly supporting the unification 
scheme. 

In this paper we focus on NGC 6251, which is 
a well studied FR I radio galaxy with a well colli- 
mated one-sided radio jet and a large scale emis- 
sion. Its total radio power is Ln^ = 1.2 x 10 32 
erg s _1 (Waggett et al. 1977), i.e. close to the 
FR I/FR II transition. However, its radio mor- 
phology is well representative of the FR I class. 
The projected linear size of the extended radio 
structure is 1.7 Mpc (for H = 75 Kms -1 Mpc -1 ), 
which makes NGC 6251 one of the largest known 
radio galaxies Perley, Bridle, & Willis (1984). The 
radio source is associated with a giant elliptical 
galaxy, which hosts a 4 — 8 x 1O 8 M black hole 
Ferrarese & Ford (1999). 

By taking advantage of the recent discovery 
that the EGRET gamma-ray source J 162 1+8203 
is associated with NGC 6251 (Mukherjee et al. 
2002), we perform a detailed analysis of its overall 
SED as we did for Centaurus A. Our result is that 
the nucleus of this low-power radio galaxy can also 
be interpreted as a misoriented BL Lac. 

The paper is organized as follows: in Sect. 2 we 
describe the HST, BeppoSAX and Chandra data 
and the analysis procedures. In Sect. 3 we deter- 
mine the SED; in Sect. 4 we discuss the nature 
of the IR-to-UV nuclear source, which is crucial 
in order to model the overall SED. In Sect. 5 we 
model the SED in the frame of a synchrotron self- 
Compton scenario, and we discuss the resulting 
model parameters, while in Sect. 6 we draw con- 
clusions. 

A redshift of z = 0.0249 for NGC 6251 (NED) 
is adopted throughout this paper. 

2. Optical and X-ray data 

In order to determine the SED, we first focus on 
the analysis of the IR-to-UV band and X-ray band 
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Fig. 1. — The central regions of NGC 6251 as they appear in HST images, from the IR to the UV band. The 
pivot wavelengths of the filters are 16060 A, 8336 A, 5461 A, 2329 A (Biretta et al. 2002). The projected 
size of the images is 10" x 10" for all images except for the UV (F220W) FOC image which it is 6" x 6". 



observations. These spectral regions are crucial in 
order to model the observed radiation, since the 
data allow us to set the position of the emission 
peaks, i.e. the frequencies at which most of the en- 
ergy is released by the source, and thus constrain 
the model. 

2.1. HST data 

We have measured the flux of the central com- 
pact source in archival HST images, from the IR 
to the UV. In Table 1 we report the details of 
the observations, and in Fig. 2 we show a picture 
gallery of the nuclear regions of NGC 6251 as they 
appear in the images. The F205W and F555W fil- 
ters might have contribution from emission lines 
which is however expected to be small (a few %), 
due to their small equivalent width compared with 
the width of the filters passband. The other filters 
are centered in spectral regions relatively free of 
strong emission lines. Therefore (with the possi- 
ble exception of F205W and F555W) the flux we 
measure is continuum emission. 

The procedure is as follows. We perform aper- 
ture photometry, evaluating the total counts in a 
circular region centered on the point source, with 
different radii for the different filters. This is im- 
portant in order to have a reliable determination 
of the total flux for the point source allowing for 
slightly different shapes of the point spread func- 
tion. For the FOC UV data both the determi- 
nation of the background and the choice of the 
aperture radius is straightforward, since the nu- 
cleus is by far the dominant component in the 
image. For the WFPC2 images, due to the pres- 



ence of the host galaxy stellar component and an 
extended dusty disk which surrounds the unre- 
solved nucleus, the largest source of error is the 
determination of the background. Thus we have 
adopted the rather well tested procedure of mea- 
suring the flux in a radius of 5 pixels (correspond- 
ing to 0.2 arcsec), measuring the background in 
the surrounding annulus (1 pixel width). This 
is essentially the same procedure as we used for 
the analysis of the complete samples of 3CR radio 
galaxies (Chiaberge et al. 1999; Chiaberge, et al. 
2002a). 

Due to the different shape of the PSF in the 
different filters, the lower resolution at IR wave- 
lengths and the larger pixel size of the detector, 
the NICMOS data have to be corrected for PSF 
loss. We have produced synthetic PSFs with the 
TINYTIM software (Krist 1995) for the different 
instrument configurations. We compared the ob- 
served shape of the nuclear PSF in NGC 6251 
with the respective synthetic PSF and we selected 
an appropriate aperture radius for the photome- 
try (r = 5.5 pixels, r = 4 pixels, r = 3 pixels 
for the F205W, F160W and F110W filters, respec- 
tively). This procedure ensures that we measure 
the background as close as possible to the nucleus 
without suffering contamination from the nucleus 
itself. However, we have determined that in this 
case the aperture photometry radius includes only 
a fraction of the total counts. Therefore, in or- 
der to estimate the total flux of the point source, 
we have derived the renormalization factors from 
the analysis of the synthetic PSFs (k — 1.234, 
k = 1.37, k = 1.38 for the F205W, F160W and 
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F110W filters, respectively) and we have corrected 
the measured flux accordingly. 

The effective frequencies of the observations 
have been calculated using SYNPHOT, by tak- 
ing into consideration the observed slope of the 
continuum emission. The count rates have been 
multiplied by the value of the header parameter 
PHOTFLAM, which represents the inverse sensi- 
tivity on the filters in the case of a flat continuum 
(F\ oc A* 3 , (3 — 0, which corresponds to F v oc v~ a , 
a = 2). We have tested that the observed slope 
does not affect the values of PHOTFLAM signifi- 
cantly (less than 1%). The fluxes measured in the 
HST observations are reported in Table 2. We es- 
timate that the typical error on the nuclear flux 
measurements is ^ 10%. The IR-to-UV slope is 
a = 1.75 ± 0.16, calculated between 16000 A and 
2200 A. 

2.2. X-ray observations and results 

We have analyzed unpublished recent BeppoSAX 
and Chandra X-ray data. These observations are 
available in the respective public archives. We fo- 
cus in particular on the BeppoSAX observations, 
since NGC 6251 is known to be a bright X-ray 
source (2-10 keV flux of ~ 10 -12 erg s _1 cm~ 2 ; 
Sambruna, et al. 1999) and therefore the Chandra 
data are significantly affected by pileup. However, 
it is extremely useful to compare the two obser- 
vations. Given its arcsec spatial resolution, which 
is ~ 200 times sharper than that of BeppoSAX, 
Chandra is indeed the ideal experiment to pin- 
point the nuclear emission. In Sect. 2.2.2 we show 
that the two observations have consistent results, 
once that the pileup effects in Chandra data are 
taken into account. 

2.2.1. BeppoSAX data 

NGC 6251 was observed in July 2001 by the 
LECS, MECS and PDS instruments on board 
BeppoSAX (with exposure times of 30.2, 80.9 and 
36.2 ksec, respectively), as part of an A05 GO 
observation (code 51322001, P.I. M. Guainazzi). 
The LECS and MECS proportional counters are 
operating in the 0.1-10 keV and 1-10 keV bands, 
respectively, while the PDS instrument is a scintil- 
lation detector operating in the high energy band 
15-300 keV. In our spectral analysis we consider 
only the range of energies in which the instruments 



are best calibrated (i.e. 0.12-4 keV for the LECS, 
1.6-10 keV for the MECS and 15-220 keV for the 
PDS). 

We perform data analysis by using standard 
procedures. We extract the source spectrum from 
an aperture of 8 and 4 arcmin radius in the LECS 
and MECS images, respectively, and we evalu- 
ate the background from blank sky fields. The 
PDS data reduction is as performed by the archive 
pipeline. We perform the spectral analysis of 
MECS and LECS data by using XSPEC vll.O on 
rebinned spectra with at least 20 counts per bin 
to use the x 2 -statistics in the fitting process. Er- 
rors are given at the 90% confidence level for one 
interesting parameter. The cross-calibration con- 
stant between LECS and MECS, obtained by fit- 
ting simultaneously the LECS and MECS data in 
the common band 1.6-4.0 keV, is found to be 0.67. 
The cross-calibration constant between the MECS 
and the PDS, constrained to be in the fiducial 
range 4 0.77-0.93, is found to be 0.88. Given the 
PDS large field of view (78 arcmin HWZI) and lack 
of spatial resolution, other sources of high energy 
photons might in principle contaminate the PDS 
data. We search for X-ray sources in the MECS 
image (25 arcmin radius FOV) and, at larger off- 
sets, in the ROSAT source catalogs 5 . Based on 
the fluxes and spectra of the X-ray sources within 
the PDS FOV, we estimate that more than 90% 
of the PDS photons are produced by NGC 6251. 

The 0.1-200 keV spectrum is well described 
(x 2 /dof = 180.5/183) by a power law with a pho- 
ton index T — 1.79 ± 0.06 mildly absorbed by a 
column density Nh = 4.7l 4 ; 7 ' x 10 20 cm~ 2 in ex- 
cess of the Galactic value (N H Gai — 5.5 x 10 20 
cm~ 2 ; Dickey & Lockman 1990). Similar results 
have also been found by Guainazzi et al. (2003). 
Since the nuclear source is embedded in a diffuse 
gas halo (e.g. Birkinshaw & Worrall 1993), we add 
in the spectral fit a thermal component using the 
mekal model of XSPEC. In this case, we obtain 
a comparably good fit (x 2 /dof = 172.6/180), al- 
though the estimate of the nuclear absorption rises 
to Nh = 1-4^0 7 x 10 22 cm~ 2 , while the power-law 
spectral index remains unchanged (r = 1. 78^34) . 

4 see the BeppoSAX handbook (Fiore, Guainazzi & Grandi, 
1999, Cookbook for BeppoSAX NFI Spectral Analysis) 
http: //www. asdc . asi . it/bepposax/sof tware/ index . html 

5 http : //www . xray . mpe . mpg . de/cgi-bin/rosat/ sre-browser 
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Table 1 
hst observations 



Instrument 


camera/config 


Filter 


HST Prop. ID 


Obs. Date 


NICMOS 


NIC2 


F205W 


G07862 


1998 Jul 06 


NICMOS 


NIC2 


F160W 


G07862 


1998 Jul 06 


NICMOS 


NIC2 


F110W 


G07862 


1998 Jul 06 


WFPC2 


PC 


F814W 


G06276 


1995 Jun 28 


WFPC2 


PC 


F555W 


G06276 


1995 Jun 28 


WFPC2 


PC 


F547M 


G06653 


1996 Sep 13 


FOC 


f/96 


F342W 


G06246 


1996 Feb 19 


FOC 


f/96 


F220W 


G06891 


1997 Jul 07 



Table 2 
HST Photometry 



Filter 


Eff. Wav. 


Flux 






A 


erg cm~ 2 s~ 


1 A- 1 


F205W 


21063 


7.3 x 10 


-17 


F160W 


16250 


6.2 x 10" 


-17 


F110W 


11777 


6.8 x 10" 


-17 


F814W 


8124 


9.1 x 10" 


-17 


F555W 


5558 


1.3 x 10" 


-16 


F547M 


5501 


1.2 x 10" 


-16 


F342W 


3435 


8.4 x 10" 


-17 


F220W 


2414 


9.1 x 10" 


-17 
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The thermal component is found to have a tem- 



perature kT = 1.20 



+0.97 
-0.34 



keV, with a metal abun- 



dance Y = 0.09l£55. The latter estimate of the 
nuclear absorption, the power-law spectral index, 
as well as the halo temperature and metal abun- 
dance, are in good agreement with the 1994 ASCA 
observation (Sambruna et al. 1999) 6 . However, 
we point out that the 2-10 keV BeppoSAX flux 
(4.78 ±0.07 x 10~ 12 erg s^ 1 cm" 2 ) is significantly 
higher than the flux measured in the ASCA obser- 
vation (1.40±0.04x 10~ 12 erg s -1 cm~ 2 ), suggest- 
ing that the source flux is variable in the X-rays. 

2.2.2. Chandra data 

NGC 6251 was observed by Chandra (37.4 ksec 
with ACIS-I) in September 2000, as part of an 
AOl GO project (Obs. ID 847, P.I. Kerp). The 
source is 3 arcmin away from the telescope aim- 
point and falls in the gap among the 4 ACIS-I 
CCDs, where the exposure map is low. This re- 
duces the number of detected photons, but does 
not mitigate the pileup effects. Based on the Bep- 
poSAX spectrum and flux, we estimated that in 
the Chandra observation (with a frametime of 3.24 
sec) the pileup fraction is above 60% (the pileup 
fraction decreases to ~ 20% if the ASCA flux is 
used). We used CIAO v2.2 to perform a standard 
data reduction and XSPEC vll.O for the spectral 
analysis. As expected, due to the high pileup frac- 
tion, the ACIS-I spectrum of NGC 6251 is much 
harder than that observed by BeppoSAX, having a 
photon index V = 1.30±0.26. By using the MARX 
4.0 simulator 7 we verified that both the spectral 
slope and the flux are in excellent agreement with 
the Chandra observation, under the assumption 
that the input spectrum is that observed by Bep- 
poSAX. We point out that it is appropriate to 
compare such observations, since the BeppoSAX 
spectrum was obtained much closer in time to the 
Chandra data than the ASCA observation. 



6 Sambruna et al 
-2 



(1999) obtained 
+0.21 
0.18 



cm"', T = 1.83^0.^8' kT = L04 ±o.i8 keV > Y = 0-06lg;^, 
while a somewhat steeper powerlaw and lower absorption 
have been found by Turner et al. (1997) in the same ASCA 
dataset 

7 Wise M.W. Davis, J.E., Huenemoerder D.P., Houck 
J.C., Dewey, D., 2003, MARX 4.0 Technical Manual: 
http : //space . mit . edu/CXC/MARX 



The nuclear Spectral Energy Distribu- 
tion of NGC 6251 
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Fig. 2. — The overall spectral energy distribution 
of NGC 6251, from the radio to the gamma-ray 
band. Details of the data are in the text. The 
dashed line corresponds to the SSC model we de- 
rive in Sect. 5. 

In Fig. 2 we show the overall nuclear SED of 
NGC 6251. The radio data are taken from Jones 
et al. (1986). The VLBA fluxes of the unresolved 
radio core are a factor of 2.5 lower than the VLA 
both at 6cm and 18cm. The IR-through-UV HST 
data and the BeppoSAX X-ray observations and 
data analysis have been discussed in Sect. 2.1 and 
Section 2.2, respectively. The EGRET spectrum 
is the fit to the spectrum of the source, as taken 
from Mukherjee et al. (2002). The SED in a 
log v — log vF v representation has two broad peaks 
and, although the peak frequencies are not well 
constrained, they are clearly located in the far IR- 
mm region and between the X-rays and gamma- 
rays. 

As in the case of Centaurus A (Chiaberge et al. 
2001), the SED appears remarkably similar to the 
typical "blazar-like" SED. In such objects, the low 
and high energy peaks are believed to be produced 
by non-thermal synchrotron and inverse-Compton 
radiation, respectively. In the following we discuss 
in detail the origin of the optical emission and we 
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show that the SSC emission model can also ac- 
count for the nuclear radiation of NGC 6251. 

4. The nature of the active nucleus in 
NGC 6251 and the origin of the opti- 
cal nuclear emission 
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Fig. 3. — The correlation between the radio core 
at 5GHz and the optical (R band) cores seen in 
HST images, for the FR I radio galaxies of the 
3CR sample (Chiaberge et al. 1999). The black 
filled circle is representative of NGC 6251. 

There is growing evidence that the nuclei of 
FR I radio galaxies are basically unobscured, and 
the compact nuclei observed in the majority of 
FR I are non-thermal radiation produced at the 
base of the relativistic jet. The main evidence for 
this is the presence of a strong linear correlation 
between FR I's optical and the radio nuclear ra- 
diation, the latter being definitely of non-thermal 
origin (Chiaberge et al. 1999; Verdoes Kleijn et al. 
2002; Capetti et al. 2002). Although the general 
behavior of the FR I class as a whole seems to 
be understood, for "single sources" a certain am- 
biguity on the nature of the optical nuclei might 
still be present. Since the IR through UV radia- 
tion is crucial for constraining the frequency of the 
synchrotron peak, it is very important to assess its 
nature. On the other hand, for the hard X-ray and 
gamma ray emission the interpretation as inverse- 



Compton emission appears to be the only viable 
possibility. 

The origin of the bright optical core of NGC 6251 
has been discussed by Ferrarese & Ford (1999). 
There are two scenarios proposed by these au- 
thors: a non-thermal process (most plausibly syn- 
chrotron emission), or thermal radiation from an 
optically thick accretion disk. 

Let us consider the position of NGC 6251 in di- 
agnostic diagrams for the nuclei of radio galaxies. 
In Fig. 3 we plot NGC 6251 in the plane formed 
by the luminosity of the optical nucleus of versus 
the radio core. The open circles are the nuclei 
of FR I radio galaxies from the 3CR catalog, as 
taken from Chiaberge et al. (1999). The dashed 
line is the regression line between the two quanti- 
ties. The black point, representative of NGC 6251, 
lies just above the correlation. This is a strong clue 
that NGC 6251 behaves exactly as other FR I. We 
stress that the comparison with the 3CR sample is 
extremely appropriate, since Laing et al. (1983), 
have included NGC 6251 in their complete revised 
3CR sample (the 3CRR). 

In Fig. 4a we show the position of the nucleus 
of NGC 6251 in the plane formed by the radio-to- 
UV vs. the optical-to-UV flux ratios. As discussed 
in Chiaberge et al. (2002b), the position of the 
objects in this plane is sensitive both to absorp- 
tion and to the different intrinsic spectral proper- 
ties of the nuclei. Steep spectrum radio quasars 
(plotted as stars) and broad lined radio galax- 
ies (empty circles) occupy the same region of the 
plane, and clearly differentiate from FR I mainly 
because of the presence of the thermal blue bump 
in the optical-to-UV spectral range. NGC 6251 
again lies in the region of FR I. 

However, as already pointed out above, for "sin- 
gle objects" we cannot exclude on the basis of the 
radio-optical properties that the optical emission 
could be scattered light produced by an accretion 
disk obscured to our line-of-sight. In this case, 
we know from Chiaberge, et al. (2002a) that its 
representative point would lie on the radio-optical 
correlation of Fig. 3 by chance. This has been in- 
deed found to be rather common among FR II in 
which strong narrow emission lines are observed. 

In order to remove the degeneracy in the in- 
terpretation of the optical emission, we use the 
[OIII]5007 emission line, which is indeed a good 
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Fig. 4. — Left panel (a): Radio core/UV nuclear flux ratio vs. optical/UV nuclear flux ratio for 3CR FR I 
(filled circles), FR II (empty circles) and radio loud quasars with z < 0.3 from the sample of Elvis et al. 
(1994). The empty circle that lies among FR I is 3C 338, a low- excitation galaxy (see discussion in Chiaberge 
et al. 2002b). Right panel (b): Equivalent width of the [OIII]5007 emission line (measured with respect to 
the optical nuclear emission) plotted vs. the ratio between the optical nuclei to radio core flux, for a sample 
of 3CR FR I and FR II (Chiaberge, et al. 2002a). Filled symbols represent detections, empty symbols are 
upper limits to the optical nuclear emission. Different symbols correspond to the different subclasses of radio 
galaxies. The shaded area represents the dispersion (ler) of the linear correlation between radio core and 
optical nuclear luminosity found for FR I. 



indicator of the intrinsic ionizing radiation field. 
From ground based observations (2.7" x 4" aper- 
ture, centered on the galaxy nucleus) the lumi- 
nosity of the [OIII]5007 line is 5 x 10 39 erg s" 1 
(Shuder & Osterbrock 1981). A useful diagnostic 
is the ratio between the flux of the emission line 
and the observed nuclear continuum emission (i.e. 
the "nuclear" equivalent width). This has been 
used to discriminate between objects in which the 
nuclear source is seen directly and objects in which 
the ionizing radiation field is obscured to our line- 
of-sight and observed only through scattered radi- 
ation (Chiaberge, et al. 2002a). In the latter case, 
high values of the equivalent width are expected, 
since only a few percent of the intrinsic nuclear 
source is observed. In Fig. 4b we show the posi- 
tion of the source in the plane formed by the [OIII] 
EW plotted against the ratio of the optical nuclear 
flux to radio core flux i? _ r . The shaded area cor- 
responds to the dispersion (la) of the correlation 
between these two quantities for FR I galaxies. 
The equivalent width of the [OIII] 5007 emission 
line for NGC 6251 is typical of the other unob- 
scured FR I nuclei. As already known from Fig 3, 



the value of i? -r is also in the same region as the 
other FR I. 

In the light of the above analysis, we conclude 
that the dominant physical origin for the optical 
nucleus is non-thermal synchrotron emission pro- 
duced at the base of the jet. Any other thermal 
contribution to the observed flux should then be 
considered lower than (or equal to) the observed 
synchrotron radiation. 

4.1. What is the role of absorption in 
NGC 6251? 

Before modeling the nuclear emission, we must 
briefly consider the possible presence of absorption 
on the line-of-sight to the nucleus, since one of the 
two models for the X-ray emission includes a cer- 
tain amount of Njj in excess of the Galactic value. 
Absorption would indeed naturally steepen the IR- 
to-UV slope of the continuum emission, and could 
in principle "artificially" produce a two-peak SED. 
This is an issue for objects like Cen A, because 
of the presence of the extended dust lane which 
certainly absorbs the nuclear radiation (Marconi 
et al. 2000). However, in NGC 6251 there is 



8 



no apparent reason for absorption in the optical. 
Firstly, the dusty disk does not appear to hide the 
nucleus. Although the extinction from the disk 
is only Ay ~ 0.6mag (in the external regions), 
this would significantly affect the UV emission 

(A ° ~ 2.5A V , Cardelli et al. 1989). If the 

v 2400A v ' ' 

disk were to cover the nucleus we would expect to 
see a progressive steepening of the spectral slope 
as frequency increases, which is not observed. Fur- 
thermore, the typical optical-to-UV spectral slope 
of FR I from the 3CR catalog lies somewhere be- 
tween a -uv = 1 — 2 (Chiaberge et al. 2002b), 
while only those objects in which absorption fea- 
tures are apparent have steeper slopes. The value 
of a .jj V = 2.0 ± 0.2 that we derive for NGC 6251 
appears to be consistent with other unobscured 
FR Is. 

In Fig. 5 we show the effects of a small 
amount of optical extinction on the broad-band 
spectral properties of NGC 6251 when compared 
with other radio galaxies. We have overplotted 
the de-reddening line for NGC 6251 on the plane 
of Fig. 4a, assuming a range of Ay from to 7. 
This confirms that only an Ay £ 1 — 1.5 mag is al- 
lowed. For higher values of extinction, the broad 
band spectral properties of NGC 6251 would be 
incompatible with any of the other sources in the 
same range of redshift and total power. In particu- 
lar, it appears that it is not possible to "drag" the 
source into the region occupied by quasars if we 
only consider the effects of obscuration. Note that 
since the typical error for each single flux measure- 
ment is £ 10%, the resulting error on the data 
plotted here is of the order of the size of the sym- 
bols. Therefore, we exclude that the nuclear emis- 
sion is thermal radiation from an accretion disk 
seen directly. 

Similarly to Fig. 5 we have overplotted the de- 
reddening curve on the EW([OIII]) vs. i? _ r diag- 
nostic plane (Fig. 6). An Ay ^1 — 2 mag would 
again move the point representative of NGC 6251 
in a region of the plane where no sources are found. 

The fact that the X-ray data can be modeled 
with absorption in excess of the galactic value, 
while in the optical there is no evidence for sig- 
nificant extinction is not surprising. In fact, the 
standard conversion between Nh and Ay, adopt- 
ing local Galactic dust to gas ratio (which in our 
case would imply Ay = 5 x 10~ 22 A ff ~ 7 mag) 
generally does not hold in AGNs. Maiolino et al. 



(2001) have found that in AGNs E(B - V)/N H 
is always lower than the Galactic value by a fac- 
tor of 3 to 100. Several different explanations for 
such a discrepancy have been proposed by various 
authors, e.g. non-standard values for the dust-to- 
gas ratio in the "absorbing" structures, a different 
location for the optical and X-ray sources and/or 
for the absorbing structures, or the presence of 
large dust grains which would flatten the extinc- 
tion curve and significantly lower the Ay /Nh ratio 
(Maiolino, Marconi, & Oliva 2001). 

As noted above, from the analysis of the X-ray 
data we cannot discriminate between the two mod- 
els that fit the data. However, for our purpose, the 
most important parameter we derive from the X- 
ray data is the slope of the (non-thermal) power- 
law component, which remains unchanged in ci- 
ther models. 

In the following we discuss the nature of the op- 
tical emission through various diagnostics. Then 
we proceed to model the overall SED in the syn- 
chrotron self-Compton scenario, with the same 
method as usually done for blazars. 

5. SSC modeling of the SED 

We model the SED of NGC 6251 in the frame 
of a homogeneous synchrotron self-Compton sce- 
nario (SSC) , which is successfully used to account 
for the overall emission of blazars (e.g. Ghisellini 
et al. 1998). Details of the model are given in Chi- 
aberge & Ghisellini (1999) and a short summary 
is reported for clarity in Appendix A. The model 
is particularly appropriate for low-power objects 
(BL Lacs), in which the external radiation field 
is probably low, and the dominant source of seed 
photons for inverse-Compton scattering is the syn- 
chrotron radiation field itself. Since the total radio 
power of NGC 6251, {L 178 = 1.2 x 10 32 erg s" 1 , 
Waggett et al. 1977) is in the range spanned by BL 
Lacs, and close to its median value (see e.g. Perl- 
man et al. 1996, Kollgaard et al. 1996, Cassaro et 
al. 1999), the use of the SSC model appears to be 
appropriate. 

In Table 3 we show the parameters of the 
model which best reproduce the SED. All of them 
are typical of BL Lacs, except for the value of 
S, the relativistic beaming factor 8 . In order to 

<5 is defined as <5 = 1/T(1 — /3cos0), where T is the bulk 
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Table 3 



Model parameters 



B = 0.03 G 

L inj = 4.0 x 10 44 erg s" 1 5 = 3.2 

7mm = 4.0 X 10 3 "iraax = 4.5 X 10 4 

P = 1.9 tesc = 15R/C 
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Fig. 5. — Radio core/UV nuclear flux ratio vs. op- 
tical/UV nuclear flux ratio for 3CR FR I (filled cir- 
cles), FR II (empty circles) and radio loud quasars 
with z < 0.3, as in Fig. 4b. The solid line is the 
dereddening curve for NGC 6251 (see text for de- 
tails). 




-6 -5 -4 -3 -2 
Log R D _ r 



Fig. 6.— Equivalent width of the [OIII]5007 emis- 
sion line plotted vs. the ratio between the optical 
nuclei to radio core flux, for the sample of 3CR 
FR I and FR II (Chiaberge, et al. 2002a), as in 
Fig. 4b. The solid line is the dereddening curve 
for NGC 6251 (see text for details). 
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perform a more detailed analysis, we compare 
our model parameters with the large BL Lac 
sample modeled by Ghisellini et al. (1998). We 
use for comparison the compactness parameter 
hnj — L in jaT/(Rm e c 3 ), which represents the 
ratio of the injected power to the size of the 
source. For NGC 6251 we obtain £ in j = 0.2, 
which lies at the higher end of the distribution 
of £inj for BL Lacs. Both the magnetic field B 
and the size of the source are well inside the range 
spanned by BL Lacs (10~ 2 £ Bbllgc ~ 1.2G and 
10 16 < R B LLac £ 10 17 cm). 

As already noted above, 5 is lower than in BL 
Lacs, where 8 ^ S £ 23 and the median value is 
16. In the analysis of the SED of Centaurus A, 
we obtained an even lower value, 5 = 1.2 (Chi- 
aberge et al. 2001). In that case, we showed that 
this could be explained by the presence of a low 
velocity component in the jet, e.g. a shear layer, 
which could be the origin of the observed nuclear 
radiation. This has been interpreted as a further 
evidence that relativistic jets have at least two ve- 
locity components: a fast "spine", with plasma 
Lorentz factors F ~ 10 — 20 which dominates the 
BL Lac emission, and a slower, mildly relativistic 
layer, which accounts for the emission of the nuclei 
of FR I radio galaxies. Considering all possible in- 
trinsic jet velocities, the maximum viewing angle 
for S = 3.2 is 9 = 18° obtained for 5 = T = 3.2. In 
NGC 6251, S = 3.2 is still compatible with "typ- 
ical" BL Lac Lorentz T factors but, in this case, 
the viewing angle to the jet axis would be even 
smaller (for T = 15 and 6 = 3.2 the viewing angle 
is 9 = 11°). 

The value for the angle of the line-of-sight to 
the jet axis we derive is rather small, but still 
compatible with the constraint 9 < 47° obtained 
from VLBI observations of the jct-counterjet ratio 
(Jones & Wehrle 2002). In absence of a firm deter- 
mination of the viewing angle to the radio jet axis, 
we cannot distinguish between the two scenarios 
proposed above for the location of the emitting re- 
gion, i.e. the fast "BL Lac" component or a slower 
(but still relativistic) region of the jet. 

Our estimate of 9 is substantially lower than 
both the orientation of the external dusty disk 
(9d,ext ~ 76°) and the inner gas disk (9 d ,inn ~ 
36°), as modeled by Ferrarese & Ford (1999). 

Lorentz factor of the jet and 6 is the viewing angle. 



These authors showed that the structure is warped 
and tilted , and the axis of the inner disk is sig- 
nificantly twisted with respect to the outer dusty 
disk. Therefore, our upper limit 9 = 18° on the 
jet direction confirms and extends their picture. 
The jet axis is not perpendicular to any of the ob- 
served disks: the warping of the accretion struc- 
ture around the central black holes would continue 
towards the innermost regions, where the inner ac- 
cretion disk and the jet are likely to be perpendic- 
ular. 

The small viewing angle 9 we find would imply 
a deprojected total size of the source ^ 5.5 Mpc, 
which would make NGC 6251 by far the largest 
known radiogalaxy. However, we point out that 
the morphology of its jets, on the largest scales, is 
complex and presents significant bending. There- 
fore, our determination of 9, which refers to the 
inner jet, is probably not appropriate to estimate 
the overall deprojected size of this radio source. 

6. Conclusions 

We have determined the overall spectral energy 
distribution of the nucleus of NGC 6251, from 
the radio to the gamma ray band, by collecting 
data from the literature and analyzing HST and 
X-ray data taken with BeppoSAX and Chandra. A 
significant improvement with respect to previous 
work is the homogeneous analysis of HST data, 
which span the range of wavelengths from the IR 
(~ 2.2^m) to the UV (~ 2500 A) and allows us 
to define a broad peak in the SED at low ener- 
gies, and the analysis of previously unpublished 
X-ray data. The X-ray data, together with the 
gamma-ray EGRET detection recently published 
by Mukherjec et al. (2002), define a further peak 
at high energies. The SED is therefore similar 
to that of BL Lac objects, in agreement with the 
AGN unification model. 

Before proceeding to model the SED, a key 
point to be addressed is the nature of the op- 
tical source. We have shown that in diagnostic 
planes for the nuclear emission of radio galax- 
ies, NGC 6251 follows the general behavior of 
FR I. Therefore, as for the large majority of 3CR 
FR I radiogalaxies, the unresolved nucleus seen in 
HST images at all wavelengths can be interpreted 
as synchrotron emission. At high energies, the 
BeppoSAX hard (a = 0.78) power-law spectrum 
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and the gamma ray radiation are most plausibly 
produced by the inverse-Compton process. 

In light of the above results, we have modeled 
the SED in the framework of a synchrotron self- 
Compton scenario, which well accounts for the 
emission of blazars in general, and especially of 
BL Lac objects. The parameters of the model are 
very similar to those found in BL Lacs, except for 
the beaming factor, as qualitatively expected by 
the unification models. More quantitatively, 5^3 
implies that the viewing angle to the jet axis can- 
not be wider than 9 ~ 20°, in agreement with the 
radio information which claims for a small viewing 
angle. 

In summary, our analysis of the overall SED 
of NGC 6251 directly supports the FR I-BL Lac 
unification scheme and strengthens the results we 
have obtained for the closest FR I, Centaurus A. 
However, it is still possible that the radiation 
which dominates BL Lacs emission is not produced 
co-spatially with the nuclear radiation observed in 
radio-galaxies. Because of the larger viewing an- 
gle to these latter sources, it is plausible that the 
"blazar" component is beamed away from the ob- 
server's line of sight, and the observed radiation 
is produced in a slower region of the jet, possibly 
a shear layer, as proposed for other sources (Chi- 
aberge et al. 2000, 2001). 

This work has been partly supported by an ESA 
post-doctoral fellowship. The authors thank the 
referee for insightful comments. MC whishes to 
thank G. Ghiscllini for useful discussion. 
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A. SSC model 



The model is described in detail in Chiaberge & Ghiscllini (1999) and Chiaberge et al. (2001), in which 
it has been applied for the first time to reproduce the nuclear SED of a radio galaxy, Centaurus A. Here we 
summarize the main assumptions of the model. 

The source is a spherical homogeneous region, embedded in a tangled magnetic field. Relativistic electrons 
are continuously injected at a rate Q(^) [cm~ 3 s _1 ] oc j~ p between 7 m ; n and 7 max (7 being the Lorentz factor), 
and they loose their energy radiatively. 

The parameters of the model are: the size of the source R, the magnetic field B, the injected luminosity 
L; n j, the relativistic beaming factor 5 = [T(l — /3cos#)] _1 (where 9 is the angle between the jet axis and 
the line of sight), 7 m i n , 7max and the slope p. The resulting electron distribution at equilibrium is a broken 
power-law. Electrons are also allowed to escape the emitting region (thus reducing the particle density, 
mimicking source expansion) in a timescale t osc . 

We obtain the equilibrium solution of the continuity equation which governs the temporal evolution of 
the emitting electron distribution N(j,t). 



where 7 = 7 s + 7c is the total (synchrotron + self Compton) cooling rate. We solve the equation numerically 
(see Chiaberge & Ghiscllini 1999) and we calculate self-consistently the produced equilibrium spectrum. We 
take into account the effects of the Klein-Nishina decline, using an approximated step function (Zdziarski 



dt 



= ^[7(7,Wy,*)] + Q(7,*)- 




= 



1986). 
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